Supernova 2007bi was a pair-instability explosion 
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Stars with initial masses 10 Mq < Mjnitiai ~ 100 Mq fuse progressively heav- 
ier elements in their centres, up to inert iron. The core then gravitation- 
ally collapses to a neutron star or a black hole, leading to an explosion — an 
iron-core-collapse supernova (SN)]!], [2]. In contrast, extremely massive stars 
(Minitiai ~ 140 Mq), if such exist, have oxygen cores which exceed M^ore = 50 Mq. 
There, high temperatures are reached at relatively low densities. Conversion 
of energetic, pressure-supporting photons into electron-positron pairs occurs 
prior to oxygen ignition, and leads to a violent contraction that triggers a 
catastrophic nuclear explosion jS], IH [5]. Tremendous energies (> 10^^ erg) are 
released, completely unbinding the star in a pair-instability SN (PISN), with 
no compact remnant. Transitional objects with 100 Mq < Mjnitiai < 140 Mq, 
which end up as iron-core-collapse supernovae following violent mass ejec- 
tions, perhaps due to short instances of the pair instability, may have been 
identified[6l [Tl H]- However, genuine PISNe, perhaps common in the early 
Universe, have not been observed to date. Here, we present our discovery 
of SN 2007bi, a luminous, slowly evolving supernova located within a dwarf 
galaxy (~ 1% the size of the Milky Way). We measure the exploding core mass 
to be likely ~ 100 Mq, in which case theory unambiguously predicts a PISN 
outcome. We show that > 3 Mq of radioactive ^^Ni were synthesized, and that 
our observations are well fit by PISN models [191 110]. A PISN explosion in the 
local Universe indicates that nearby dwarf galaxies probably host extremely 
massive stars, above the apparent Galactic limit [llj . perhaps resulting from 
star formation processes similar to those that created the first stars in the 
Universe. 



We discovered a new optical transient (SNF20070406-008) on 2007 April 6.5 (UT 
dates are used throughout this paper) at right ascension a = 13''19'"20.2* and declination 
6 = 08°55'44.0" (J2000). Follow-up spectroscopic observations showed that this was a 
supernova (SN 2007bi), with no trace of either hydrogen or helium, leading to a Type Ic 
classification [12], albeit of a peculiar nature with only one previously known counterpart 
(SN 1999as[13]; Fig. 1). No signs of interaction with circumstellar material (CSM, a 
major source of uncertainty in the analysis of previous luminous SNe[lll [7]) are observed 
throughout the evolution of this event (Fig. 1). A search for pre-explosion data recovered 
observations by the Catalina Sky Survey (CSS[T5]) that allowed an accurate determination 
of the date and magnitude of the supernova brightness peak. Photometric observations 
continued for 555 days during our intense follow-up campaign. A red (i?-band) light curve 
is plotted in Fig. 2 (top). 

The measured light curve is unique, showing a very long rise time to peak (~ 70 days; 
Fig. 2; SI § 2), an extreme luminosity reaching an absolute peak i?-band magnitude 
of Mr = —21.3 mag, and a slow decline (0.01 mag day~^ over > 500 days), consistent 
with the decay rate of radioactive ^^Co. These properties suggest that the very massive 
ejecta were energized by a large amount of radioactive nickel (> SM©; Fig. 2, 3; SI § 3), 
as expected from pair-instability SN models|H \5\ HU]. Our spectra, lacking any signs of 
hydrogen or helium, indicate that this mass is dominated by C, O, and heavier elements. 
The large amount of kinetic energy released, -Ek ~ 10^^ erg (Fig. 2; SI § 3), is comparable 
to those derived for the most energetic gamma-ray bursts (GRBs[Tn]), placing this event 
among the most extreme explosions known. In Fig. 2 (bottom) we show theoretical light 
curves calculated from PISN models[5], [9] prior to our discovery. The data fit the models 
very well, suggesting we observed the explosion of a star with a helium core mass around 
100 Mq. 

PISN models imply that such an explosion would synthesize 3-10 M© of radioactive 



56]\jj (^Xable 1). Such a large amount of newly synthesized radioactive material would 
energize the SN debris for an extended period of time, ionizing the expanding gas cloud. 
Collisional excitation would lead to strong nebular emission lines, whose strength should 
be roughly proportional to the amount of radioactive source material, providing another 
testable prediction. Fig. 3 (top) shows a comparison of the nebular spectrum of SN 
2007bi with that of the well-studied, ^^^Ni-rich SN 1998bw, which produced ~ 0.5 Mq of 
radioactive nickel [T7]. suggesting that SN 2007bi produced ^ TM© of nickel (SI § 3), again 
supporting the PISN interpretation. 

Modelling the nebular spectrum, we are able to resolve the elemental composition of 
the fraction of the ejected mass that is illuminated by radioactive nickel. We can directly 
measure the abundances of C, O, Na, Mg, Ca, and Fe, and derive the mass of radioactive 
^^Ni. Our elemental abundance ratios are in good agreement with model predictions [5J for 
heavier elements, while lighter elements (C, O, Mg) seem to be underobserved. Adopting 
the calculated model output [5] for elements which do not have strong nebular emission in 
the optical (mostly Si and S, and some Ne and Ar), we arrive at a total illuminated mass 
of > 50 M0, with a composition as described in Table 1. Note that this falls well below the 
total mass derived from the photometry, indicating that even the unprecedented amount 
of radioactive nickel produced by SN 2007bi was not sufficient to energize the entire mass 
ejected by this extreme explosion (SI § 6). The unilluminated mass probably contains 
more light elements that originated in the outer envelopes of the exploding star, and seem 
deficient in our nebular observations (see SI § 3,6 for additional details). 

Our data thus provide strong evidence that we have observed the explosion of a helium 
core with M ^ 100 Mq, which, according to theory, can only result in a PISN[3l HI El O 
[inillS]- The measured light curve, radioactive nickel yield, and elemental composition of 
the ejecta are consistent with models of PISNe that were calculated before our discovery. 
Based on fewer observations of SN 2007bi, combined with their analysis of the host-galaxy 



properties, ref. [T^] consider both a PISN model and a massive iron-core-collapse SN 
interpretation [2D], slightly favoring the latter. However, our quantitative estimate of the 
helium core mass from our peak light-curve shape and the analysis of the nebular spectra, 
is inconsistent with iron-core-collapse models[20] and theoretically requires a PISN[5| [TO]. 
We thus conclude that we have most likely discovered the first clear example of a PISN. 

There are several implications of this discovery. Theory allows stars as massive as 
1000 M0 to have formed in the very early Universe [21]. However, the most massive stars 
known in the local Universe (e.g., luminous blue variables) have estimated masses around 
ISOMqIH]. In the single example known so far, such a hypergiant star exploded in a 
regular core-collapse event (SN 2005gl)[2j. Our detection of a PISN from a ~ 100 Mq 
core suggests a progenitor with an estimated initial mass around 200MQ[5j, assuming 
very low mass-loss rates appropriate for zero-metallicity stars. We note that this estimate 
is highly sensitive to poorly understood mass loss, and that high-metallicity mass-loss 
prescriptions would require an even higher initial mass. In a sense, our discovery of such 
a core is in conflict with commonly used mass-loss calculations, which do not allow such 
a high- mass core to form at the measured metallicities[19]. Our finding probably requires 
the modification of mass-loss paradigms, perhaps through increased clumping in massive 
star winds[TBl El E], at least during the hydrogen-rich mass-loss phase. 

Regardless of the exact mass loss adopted, our data indicate that extremely massive 
stars above the Galactic limit (M > 150 Mq) are formed in the local Universe. Perhaps the 
dwarf, metal-poor host galaxy of SN 2007bi {Mb = —16.3 mag at 2 = 0.1279, indicating 
an approximate metallicity of 12+log[0/H] = 8.25 [ref. [22]]; see ref. [19j for additional 
details) retained conditions that were similar to those prevalent in the early Universe. 
Luminous events like SN 2007bi can therefore serve as beacons, focusing our attention 
on otherwise unremarkable local dwarf galaxies that can be used as fossil laboratories to 
study the early Universe. 



Our observational confirmation of PISN models supports their use to predict the de- 
tectability and observed properties of PISNe from the first stars by future missions such 
as the James Webb Space Telescope, to estimate their contribution to the chemical evolu- 
tion of the Universe[5], and to calculate their impact on the re-ionization of the Universe. 
With the advent of new wide- field surveys such as the Palomar Transient Factory [231 121] 
and CRTS [15] that monitor millions of nearby low- luminosity, anonymous galaxies [25], 
many additional such events should soon be discovered and will further illuminate these 
important questions. 
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Element [Mq] 


C 


O 


Ne 


Na 


Mg 


Si 


S 


Ar 


Ca 


56Ni 


Total 


Measured (1) 


1.0* 


10.1* 


4.0 


0.0012* 


0.068* 


22.0 


10.0 


1.3 


0.75* 


4.5* 


53.8 


Measured (2) 


1.0* 


12.0* 


4.0 


0.0013* 


0.095* 


22.0 


10.0 


1.3 


0.90* 


5.7* 


57.1 


Measured (3) 


1.2* 


14.6* 


4.0 


0.0018* 


0.13* 


22.0 


10.0 


1.3 


1.00* 


7.4* 


61.7 


Measured (4) 


1.0* 


7.5* 


4.0 


0.0015* 


0.065* 


22.0 


10.0 


1.3 


0.95* 


3.65* 


50.6 


Measured (5) 


1.0* 


9.1* 


4.0 


0.0018* 


0.085* 


22.0 


10.0 


1.3 


1.10* 


4.6* 


53.3 


Measured (6) 


1.0* 


11.3* 
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0.0023* 


0.12* 


22.0 


10.0 


1.3 


1.00* 


6.0* 


56.8 


95 M0 model|5] 


4.1 


45.2 
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0.0028 
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21.3 
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100 Mq model[5] 
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43.9 


4.1 
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23.1 
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5.82 
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105 Mq modeim 
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42.7 
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4.40 


24.5 
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1.40 


9.55 
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95 M0 modelfTO] 


0.28 


38.1 
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2.32 
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100 Mq model [TOl 
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1.51 
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2.45 


22.6 


15.0 


2.54 


2.23 


8.00 


100 


105 M0 model[IO] 


0.23 


35.2 


1.47 


0.0009 


2.47 


23.9 


15.9 


2.73 


2.46 


11.89 


105 



Table 1: Predicted PISN ejecta composition compared with our measurements. We report 
six different estimates based on the two available late-time spectra (Fig. 3): VLT day 414 
([l]-[3]), and Keck day 530 ([4]- [6]). Each spectrum is modeled assuming three different 
explosion dates: 45 days ([1], [4]), 77 days ([2], [5]), and 113 days ([3], [6]) prior to peak 
(observed). The results are qualitatively similar, with ^^Ni masses of 3.7-7.4 Mq and total 
masses 51-62 Mq. Elements noted by an * mark are directly constrained from optical 
nebular emission lines. Abundances of other elements are probed through their cooling 
effects via lines outside of the optical range (which are modelled), but we consider these 
constraints to be weaker. 



Figure 1: 

Spectra of the unusual Type Ic SN 2007bi. We observed SN 2007bi on 2007 April 15.6 
and 16.4, using LRIS[26J mounted on the Keck I 10m telescope. Narrow emission lines 
(see Fig. 3 for details) indicate a redshift of z = 0.1279. A survey of our databases of SN 
spectra shows that this event is similar only to a single previous example, SN 1999as[T3]. 
until now the most luminous known SN Ic by a wide margin. SN 2007bi has a comparable 
luminosity (Fig. 2). We identify the most prominent features (marked) as arising from 
calcium, magnesium, and iron, and derive a photospheric velocity of 12,000 km s~^. 
A model fit (SI § 5) confirms these line identifications and shows that the absorption 
in the blue-ultraviolet part of the spectrum is dominated by blends of iron, cobalt, and 
nickel lines. The shallow, poorly defined trough seen around 5500-6000 A could arise from 
blends of silicon and sulfur lines, while lines of neutral oxygen and sodium that are usually 
prominent in SN Ic spectra are remarkably weak. No hydrogen or helium lines are seen, 
confirming the Type Ic classification, and strongly disfavoring the possibility that CSM 
interaction contributes significantly to the remarkable luminosity |14[ [TJ No narrow 
sodium absorption is seen, indicating little absorption by dust in the host galaxy [19]. 
The strong emission line around 7300 A seems to arise from [Ca II] emission at zero 
expansion velocity, and is usually observed in SNe only in late-time nebular spectra. A 
complete analysis of our full set of photospheric spectra will be presented in a forthcoming 
publication; see also additional details in ref. [T9] . 




Figure 1: 



Figure 2: 

Radioactive ^^Ni and total ejected mass from the light-curve evolution of SN 2007bi are 
well fit by PISN models, (a) The i?-band light curve of SN 2007bi. Data with standard 
deviation errors are shown. We have compiled observations obtained by the 48-inch (1.2 m) 
Oschin Schmidt, the 60-inch (1.5 m) robotic, and the 200-inch (5 m) Hale telescopes at 
Palomar Observatory, as well as photometry from the Catalina Sky Survey (CSS [15]) and 
synthetic photometry integrated from our late-time Keck spectrum (Fig. 3; see SI § 1,2 
for additional details). We find a peak magnitude of M/j = —21.3 ± 0.1 mag on 2007 Feb. 
21 (SI § 2). Our error is dominated by the absolute zero-point calibration uncertainty. 
The outstanding peak luminosity of this event, if radioactively driven, suggests that a 
remarkable amount of ^®Ni was produced (> 3Mq; ref. [28]; SI § 3). The slow rise 
time derived from our fit (77 days; SI § 2), combined with the measured photospheric 
velocity (12,000 kms~^; Fig. 1), requires very massive ejecta (Mej ~ 100 Mq) and a 
huge kinetic energy release [E]^ ^ 10^'^ erg; SI § 3), adopting the commonly used scaling 
relations [251 [25] . An independent direct estimate for the ^^Ni yield is obtained from 
the luminosity during the late-time decay phase, compared to that of SN 1987A[30J (SI 
§ 3) . Given the uncertainty in the explosion date of SN 2007bi and a range of bolometric 
correction values (SI § 2), the ^^Ni mass produced by SN 2007bi was 4Mq < Mse^i < 
7M0. The total radiated energy we measure by direct integration of the light curve is 
Erad ~ (1 ~2) X 10^^ erg (SI § 3), comparable to that of the most luminous SNe known]?], 
(b) Comparison of the observations of SN 2007bi with models calculated before the SN 
discovery [5l [9]. The curves presented are for various helium cores (masses as indicated) 
exploding as PISNe, and cover the photospheric phase. The data are well fit by 100- 
110 Mq models. At later times, the emission is nebular and bolometric corrections used 
to calculate the model i?-band light curve cease to hold (SI § 4). In comparing with these 
restframe models, cosmological time dilation for z = 0.1279 has been taken into account. 
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Figure 3: 

Ejecta composition from nebular spectra of SN 2007bi. (a) Shown are two late-time 
spectra of SN 2007bi. A spectrum obtained with the F0RS2 spectrograph mounted on the 
ESO 8.1m Very Large Telescope at Paranal Observatory on 2008 April 10 (414 days post- 
peak; 367 days restframe) is not completely nebular yet. Prominent broad SN emission 
peaks from neutral and singly ionized elements are marked. Multiple narrow host-galaxy 
emission lines ([O II] A3727, H/3, [O III] AA4959, 5007, Ha, [S II] AA6716, 6731) are 
seen at z = 0.1279. An even later spectrum (530 days post-peak; 470 days restframe) 
was obtained with LRIS mounted on the 10 m Keck I telescope on 2008 August 4. This 
spectrum is fully nebular, though of lower signal-to-noise ratio. Comparison with a late- 
time spectrum of the ^^Ni-rich SN 1998bw, which produced 0.5 of ^®Ni[T7], adjusted 
for the larger distance and later restframe spectroscopic observations of SN 2007bi and 
multiplied by eight, provides a good fit for intermediate-mass elements (O, Na, Mg, Ca) 
but underpredicts the strength of the Fe lines. Assuming that emission-line luminosity 
scales with the mass of energizing ^^Ni, we derive from the scaling factor a ^^Ni mass 
7.7 Mq < MseNi < 11.3 Mq (SI § 3), in reasonable agreement with estimates from early 
(peak) and late-time (radioactive tail) photometric estimates. We note that the lack of Ha 
emission at these late epochs is an especially strong argument against CSM interaction. 
A lower limit of ~ 5 x 10^^ cm on the distance to any H-rich material (in particular, 
recent mass loss) is derived from this non-detection assuming an expansion velocity of 
12,000 kms-i (Fig. 1). 

(b) Modelling the nebular spectrum. Employing our nebular spectroscopy code^/Tj we are 
able to constrain the composition of the ejecta. We find that the emission-line luminosity 
requires an initial ^^Ni mass of 3.7-7.4 Mq and the composition given in Table 1. As 
can be seen there, our measurements are well fit by the theoretical predictions of PISN 
models[5] calculated well before we discovered SN 2007bi. 
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Supplementary Information 



(1) Technical observational details 

Photometry: 

Discovery and follow-up observations of SN 2007bi were obtained using the Palomar- 
QUEST camera mounted on the 48" Oschin Schmidt telescope at Palomar Observatory 
(P48) as part of the SN Factory (SNF) program |31]. These i?-band observations were 
pipeline-reduced by the SNF software, including bias removal, fiatfield corrections and an 
astrometric solution. 

Observations using the robotic 60-inch telescope at Palomar Observatory (P60) were 
pipeline-processed [32], including trimming, bias removal, fiatfield corrections, and an as- 
trometric solution. 

Observations using the 200-inch Hale telescope at Palomar Observatory (P200) were 
obtained using the Large Format Camera (LFC) in SDSS r-band and were cross-calibrated 
onto standard R band as detailed below. 

Observations by the Catalina Sky Survey (CSS [15]) were obtained using the CCD cam- 
era mounted on the 0.7 m Catalina telescope. These unfiltered data were cross-calibrated 
onto the i?-band grid as detailed below. 

Synthetic photometry was derived from the late-time Keck spectrum using the meth- 
ods of ref. [33] . 

Table 3 provides the full list of photometric data. 

Spectroscopy: 



Early-time spectroscopy presented in Fig. 1 was obtained with the Low Resolution 
Imaging Spectrometer (LRIS[26]) mounted on the 10m Keck I telescope on Mauna Kea, 
Hawaii. The presented spectrum was obtained on Apr. 16, 2007 in long-slit mode. The 
exposure time was 600 s at airmass 1.08 under clear sky conditions with variable seeing 
around 2". The D560 dichroic was used with the 600 line mm~^ grism on the blue side, and 
the 400 line mm~^ grating blazed at 8500 A on the red side, with the 1.5" slit oriented at 
the parallactic angle [53]. The spectral resolution achieved was ~ 9A on the red side and 
~ 5 A on the blue side. Shutter problems cast doubt on the reliability of the absolute flux 
calibration, so we have forced the spectral shape to match that of a lower signal-to-noise 
ratio spectrum obtained on the previous night (April 15, 2007) using the same instrument 
in spectropolarimetry mode and an otherwise identical setup, for which a reliable flux 
calibration was obtained. 

The spectrum of SN 1999as shown in Fig. 1 was obtained ~ 3 weeks after discovery. 
Details about these data will be presented elsewhere, see initial report in ref. 

Late-time LRIS spectra presented in Fig. 3 were processed using a pipeline developed 
by one of us (MS [37]). Following standard pre-processing (e.g., overscan subtraction), the 
data are divided by a normalized fiatfield removing pixel-to-pixel sensitivity variations 
and correcting for the different gains of the CCD amplifiers. Cosmic rays are identified 
and removed using LACOSMIC[38j. We perform sky subtraction by subtracting a two- 
dimensional sky frame constructed from sub-pixel sampling of the background spectrum 
and a knowledge of the wavelength distortions as determined from two-dimensional com- 
parison lamp frames. We also perform a fringe-frame correction on the red side. The 
two-dimensional frames are transformed to a constant dispersion using comparison lamp 
exposures. The spectral extraction is performed by tracing the object position on the 
CCD and using a variance-weighted extraction in a seeing-matched aperture. An error 



spectrum from the statistics of the photon noise is also extracted. The wavelength cal- 
ibration of each extracted spectrum is then adjusted slightly using the position of the 
night-sky lines to account for any drift in the wavelength solution. Telluric correction 
and relative flux calibration is performed using spectrophotometric standard stars. We 
then finally combine the two sides into a single calibrated spectrum. We match the flux 
across the dichroic by defining narrow box filters on either side of the dichroic, and then 
use a weighted mean to combine the spectra, in the process re-binning to a constant 2 A 
per pixel resolution. The result is a contiguous object spectrum, together with an error 
spectrum representing the statistical uncertainties in the flux in each binned pixel. 

The late-time spectrum from the Very Large Telescope (VLT) presented in Fig. 3 was 
taken on April 10, 2008 using the F0RS2 spectrograph with the 300V+20 grism. Four 
3600 s exposures were obtained with the slit oriented along the parallactic angle. The 
spectra were pre- reduced (bias and fiatfield corrected), extracted, and wavelength- and 
fiux-calibrated using standard tasks within IRApy. The wavelength and flux calibration 
was computed using comparison lamps and the standard star Feige 56 observed with the 
same instrumental configuration. Telluric features were removed using the standard-star 
spectrum (observed at similar airmass). The combined spectrum revealed contamination 
from the host galaxy that was subsequently removed using an extracted spectrum on the 
edge of the galaxy where the SN contribution was negligible. 

The absolute flux level of the late VLT spectrum was calibrated with photometry 
obtained during the same night [19]. The late-time Keck spectrum does not have con- 
temporaneous photometry. However, the fact that the derived synthetic photometry is 
consistent with a smooth extrapolation from the previous photometry, and that the neb- 

-'^IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy 
Observatories, which are operated by AURA, Inc., under cooperative agreement with the National Science 
Foundation 



ular spectral analysis of both the photometry-calibrated VLT spectrum and the later 
Keck spectrum give very similar results, suggests that the absolute flux level of the Keck 
spectrum is reasonable. 



(2) SN 2007bi photometry 



We placed our i?-band P60 photometry onto an absolute grid using SDSS photometry 
of multiple nearby sources, and solving for the zero-point offset and color corrections 
for individual images using a least-squares-based solver (E. O. Ofek, in preparation). A 
similar method was employed for the P200 r-band images. Photometry of the P48 data 
was obtained using the SNF pipeline and the derived -R-band magnitude agrees well with 
other data. CSS photometry was calibrated onto an _R-band grid by applying a constant 
zero-point offset, and overlapping data agree well with Palomar photometry (Fig. 1). 
Comparison with standard-star photometry from ref. [19] shows an offset of ~ 0.1 mag 
which we attribute to residual differences due to telescope and filter transmissions. Table 3 
provides our full set of i?-band data. Additional photometry will be given in a forthcoming 
publication. 

In view of the lack of evidence for extinction of SN 2007bi in its dwarf host [19], 
including non-detection of Na D absorption lines in our spectra, we estimate that host 
extinction is negligible and correct only for Galactic extinction of A^i = 0.07mag[39] 
taken from NED. We use a distance modulus appropriate for z = 0.1279 in the standard 
cosmology [Hq = 71 kms^^ Mpc^^; = 0.27; Qa = 0.73) of m - M = 38.86 mag. 

We subtract the host-galaxy contamination from our photometry in the following 
manner. We calculate the i?-band host-galaxy magnitude i?gai = 22.54 mag using available 
SDSS measurements (rgai = 22.62 mag, Qg^i = 22.45 mag) from which we derive (using the 
formulae of ref. [^) V — R = —0.0188 mag, r — R = 0.083 mag, and therefore i?gai = 
22.62 — 0.083 = 22.54 mag. We then subtract the host-galaxy flux from all photometric 
data. 

Since the redshift of SN 2007bi is non-negligible, one needs in principle to correct the 
observed i?-band photometry into restframe /2-band photometry (so-called i^-correction. 



Kjiji). However, in view of the fact that ref. [19] shows Kji^ < 0.16 mag at all epochs, 
we neglect this correction. 

To compare with models and other SNe it is useful to convert observed filtered pho- 
tometry into bolometric photometry. Combining our optical data, early-time IR data 
(to be presented elsewhere) and data from ref. [19], we find that the bolometric correc- 
tion {BCr — Mboi — Mr) is typically ~ —0.5 mag and is always > —0.75 mag. When 
appropriate, we use the range > BCr > —0.75 mag in our calculations. 

To estimate the SN peak magnitude and rise time we fit our sparse carly-timc pho- 
tometry with low-order polynomials. Experimenting with various polynomials and data 
ranges used in the fits, we derive typical values for the rise time (defined here as the 
time required for the SN to rise by 5 mag to peak) of t^ise — 77 days and a peak mag- 
nitude Mr = —21.3 ± 0.15 mag. We show typical fits in Fig. 4 (red and blue curves), 
as well as the fit that yields our shortest rise time (45 days). Of course, we cannot 
constrain more complex light curve shapes given the sparsity of the data. We adopt 
^risc = 77 days as our fiducial rise-time value, and use the range 45 < trise < 110 days 
(restframe 40 < t^se < 97.5 days) when appropriate (110 days is a typical value derived 
from PISN models. Fig. 2b). 
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Figure 4: 



Polynomial fits to the early light-curve data. Polynomials are fit to the data in a 
least-squares sense. Due to the paucity of the data, higher order (> 5) polynomials are 
underconstrained. 



(3) Derived physical properties 



We derive the physical properties of the explosion from the observations using several 
independent methods. Our results are summarized in Table 2. 

^^Ni mass: 

We estimate the synthesized ^^Ni mass from the observed peak magnitude using the 
peak luminosity vs. ^^Ni mass correlation [2 8] and find MseNi — 3.5 M©. 

An independent estimate for the ^^Ni mass is obtained from the luminosity of the 
radioactive tail of the light curve of SN 2007bi, which closely follows the theoretical decay 
slope (0.0098 mag day"^) of ^^Co (the radioactive daughter product of ^^Ni). We used 
a direct comparison with the measured bolometric light curve of SN 1987A which is 
known to have been powered by the decay of 0.07 Mq of ^^Co at these stages[30] (Fig. 
5). We derive Msbn; = 5.3 M©, with a range of 4.4 < Mssn; < 7M0 corresponding 
to restframe rise times 40 days < tnse < 97.5 days and bolometric corrections —0.75 mag 
< BCr < Omag (see above). 

Modelling of the nebular spectra of SN 2007bi (Fig. 3b, see also below) yields estimates 
of the initial ^^Ni mass since continued radioactivity energizes the expanding ejecta and 
excites the observed emission lines. Models of the two available late-time spectra assuming 
a range of explosion dates give estimates of 3.7 Mq < Mse-^i < 7.4 Mq (see main text, Table 

!)• 

Scaling our latest spectrum to nebular spectra of the well-studied SN 1998bw (Fig. 
2 a) provides an estimate of the ^^Ni mass assuming that the strength of the emission 
lines from the intermediate-mass elements (O, Mg, and Ca) is proportional to MseNj. 
Correcting for the different restframe epochs of the spectra - times since explosion of 



SN 1998bw (373 days) and SN 2007bi (470 days) - we get an estimated ^^Ni mass range 
of 7M0 < MssNi < 11 Mq corresponding to restframe rise times in the range 40 days 
< trise < 97.5 days and assuming that SN 1998bw produced 0.5 Mq of 56Ni[17]. Note that 
this scahng is rough given the different 0/Fe ratios of SNe 1998bw and 2007bi. 

Inspecting theoretical models [5, 9] (Fig. 2b), we find that our measured light curve 
is best fit by models producing 3 M0 < Mse^i < 11 M0 (total helium core masses 
95 M0 < Mne < 110 Mq). 

Total ejected mass: 

Modelling of the nebular spectra of SN 2007bi (Fig. 3b, see also below) yields estimates 
of the total ejected mass of 51 M0 < Mgj < 61 M0 (see main text. Table 1). Note that 
these are strictly lower limits on the total mass for the following reasons. The velocity 
of the material emitting the nebular emission lines (fneb ~ 5, 600 km s~^) is much lower 
than the velocities observed during the very extended (> 100 days) photospheric phase 
(vph = 12,000kms~^; Fig. 1, see also Fig. 16 of ref. [19]. This high-velocity material 
is therefore not contributing to the spectrum at late times, and represents a substantial 
additional reservoir of mass in addition to our derived values. In addition, almost all 
models of massive He cores predict that the outer few solar masses of material will be 
composed mostly of He, while we have not observed He lines at any epoch. This is usually 
explained by the fact that He lines are non-thermally excited, and would only appear when 
^^Ni is mixed into the He envelope providing a hard non-thermal exciting spectrum |41]. 
The lack of observed helium in SN 2007bi is therefore consistent with a spatial segregation 
between an external He envelope and an internal core of heavier elements into which newly 
synthesized ^^Ni is mixed. The mass of this outer envelope (which should be several solar 
masses) should again be added to the spectroscopic estimate given here. See below (§ 6) 



for additional discussion of the ejecta geometry. 

We can derive the total ejected mass Mgj in the explosion from commonly used[28,29] 
scaling relations. The required measurements are the rise time t^se and the photospheric 
velocity fph- We note that these are rough estimates that depend on the object used 
to anchor the scaling. Here, we compare SN 2007bi with two well-studied SNe Ib/c for 
which the rise time is known (owing to their coincidence with GRB/X-ray- flash events 
that fix the explosion time): SN 1998bw and SN 2008D. We adopt Wph = 20,000kms-\ 
trise = 17days, and Mej = IIM0 for SN 1998bw[17]; Vp^, = 10,000kms-^ t^ise = 19days, 
and Mej = 7Mq for SN 2008D[42J; and Wph = 12,000kms-^ (Fig. 1) and trise = 66 
restframe days (with a range 40 days < t^se < 97.5 days) for SN 2007bi (see above). Using 
the relations from ref. [29] we get an estimated mass of Mgj ~ 105 Mq (with a range 
36 Mq < Mej < 173 Mq for restframe rise times 40 days < t^isc < 97.5 days) and almost 
exactly the same results using either SN 1998bw or SN 2008D. We note that the lower 
range of these values is inconsistent with the spectroscopic estimates provided above, 
arguing against such shorter rise times (i.e., indicating that trise > 60 days). 

Inspecting theoretical models[5,9] (Fig. 2b), we find that our measured light curve is 
best fit by models with total helium core masses 95 < Mhg < 110 Mq. 

Kinetic energy: 

We can derive the kinetic energy E]^ generated by the explosion from commonly 
used [28, 29] scaling relations. The required measurements are the rise time tnse and the 
photospheric velocity Vph- We note that these are rough estimates that depend on the 
object used to anchor the scaling. Here, we compare SN 2007bi with two well-studied 
SNe Ib/c for which the rise time is known, as discussed above: SN 1998bw and SN 
2008D. We adopt Vph = 20,000kms-\ tnsc = 17 days, and = 30 x 10^^ erg for SN 



1998bw[17]; v^^. = 10,000kms-\ t^se = 19 days, and = 6 x lO^^erg for SN 2008D|12]; 
and fph = 12, 000 km (Fig. 1) and trise = 66 restframe days (with a range 40 days 
< trise < 97.5 days) for SN 2007bi (see above). Using the relations from ref. [29] we get an 
estimated energy of E]^ ^ 115 x 10^^ erg (with a range 36 x 10^^ erg < E]^ < 273 x 10^^ erg 
for restframe rise times 40 days < tiisc < 97.5 days) and again, similar results using either 
SN 1998bw or SN 2008D. We note that the upper range of these values appears to exceed 
the total budget of available nuclear energy in a PISN (i? < 80 x 10^^ erg) [5], arguing for 
shorter rise times. It therefore appears that our fiducial rise time (observed t^se ~ 77 days; 
trise = 66 days at restframe) provides a reasonable fit considering all available constraints. 

As a sanity check we can directly estimate the kinetic energy using E]^ = (l/2)Mej v"^, 
where v is the mass-averaged expansion velocity. Assuming our fiducial estimated mass 
Mgj = 100 Mq, and that ~ 1/2 of that mass lies at low (nebular) velocities based on our 
nebular analysis while the rest travels at velocities around the photospheric values, we 
can adopt as the mean v = 8, 000 km s~^ and derive E^^ ^ 80 x 10^^ erg. 

Radiated energy: 

Direct integration under the observed light curve provides an estimate of the total 
radiated energy of 1 x 10^^ erg < i^^ad < 2 x 10^^ erg for a range of bolometric corrections 
—0.75 mag < BCn < Omag (see above). This is comparable to the total luminosity of 
the brightest SNe known[7] (see immediately below). 



Quantity 


Method 


Value [range] 


Assumptions 


^^Ni mass 


Peak magnitude 
SN 1987A comparison 

Nebular modelling 
SN 1998bw comparison 

Light-curve models 


3.5 Mq 

5.3 [4.4 .. 7] M0 
[3.7 .. 7.4] Mq 
8.9 [7.7 .. 11.3] M0 
[2.7.. 11] M0 


trise = [45. .110] days, BCr = [-0.75.. 1] mag 
trise = [45. .110] days 
trise = [45.. 110] days 
ref. [5] 


Ejected mass 


Nebular modelling 
Light-curve scaling 
Light-curve models 


> 50 Me 
105 [37 .. 173] Me 
[95 .. 110] Me 


trise = [45.. 110] days 
trise = [45.. 110] days 
ref. [5] 


Kinetic energy 


Light curve scaling 
(l/2)Mej X t;2 


132 [68 .. 273] lO^i erg 
80 10^1 erg 


trise = [45.. 110] days 
Mej = WOMq,v = 8,000kms-i 


Radiated energy 


Direct integration 


[1 .. 2] IQSI erg 


BCr = [-0.75.. 1] mag 



Table 2: Summary of physical properties derived 



(4) Comparison with other luminous SNe 

Several very luminous SNe were reported in the last few years, and have been specu- 
lated to be PISNe. Most prominent were SN 2006gy[l3l El |15] (Fig. 5), SN 2005ap|l6], 
SN 2006tf[47], and SN 2008es[27],[lH]. Two major differences distinguish SN 2007bi from 
these previous cases and strongly suggest it was a PISN. First, in all previous cases the 
spectra of these luminous SNe showed evidence for hydrogen, which is lacking in SN 
2007bi. This has a fundamental implication, since in the case of SN 2007bi, all of the 
observed ejecta had to come from the helium core, allowing us to directly constrain its 
mass, and ultimately to provide compelling evidence for a PISN (a helium core mass 
above 50 Mq). In contrast, the ejecta mass in other SNe may be dominated by hydrogen, 
complicating an attempt to constrain the helium core mass. 

Next, some of the previously studied luminous SNe showed strong signatures of CSM 
interaction Hg, 07], [7], [14], [27], which SN 2007bi lacks. Thus, the luminosity of SN 2007bi 
reflects directly on the physics of the explosion (radioactive element synthesis, explosion 
energy). The luminosity in other cases may be dominated by conversion of kinetic en- 
ergy from a more standard SN explosion into luminosity via shocks launched following a 
collision with a massive CSM, previously lost from the progenitor star. Thus, while the 
accumulated data for SN 2006gy seem to converge on an extremely massive progenitor 
(M = 100 M0 or more) [7,8], they disfavor a PISN as the underlying energy source. In 
other cases (SN 2005ap|in]; SN 2008es[lH],[27]), signatures of interaction are less clear, 
but CSM interaction is still favored as the source of the observed luminosity [2 7]. 

Compared to other hydrogen- and helium-deficient SNe of Type Ic, SN 2007bi is 
by far the most luminous and energetic, with the exception of SN 1999as[13].|5B].[lU] 
which appears quite similar, but lacks observations (especially late-time photometry and 
spectroscopy) of similar quality. Even the most energetic "hypernovae" associated with 



GRBs[5ni ED IS2] pale in comparison, with approximately an order of magnitude less 
energy released and radioactive ^^Ni produced[17] (Fig. 5). Indeed, the term "hypernovae" 
seems better suited to the truly extreme PISN explosions we have identified here [53]. 
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Figure 5: 

Comparison with other luminous SNe. Wc compare the i?-band hght curve of SN 2007bi (red circles) with that of other 
luminous SNe observed at late times: the very luminous SN 2006gy[8,14] (blue stars) and the prototypical GRB/hypernova 
SN 1998bw| 50l 1541 1551 1561 (magenta squares). We anchor our discussion on the well-studied SN 1987A whose bolometric 
light curve[30] is also presented (black dots). SN 1987A produced ~ 0.07 Mq of ^''Ni. The decay of its daughter nucleus 
^^Co drives the late-time emission, as can be seen from the comparison with the theoretical decay line plotted as the dashed 
black curve. While the relatively massive and slow ejecta of SN 1987A generally provide an efficient envelope to trap the 
radioactive energy and convert it into radiation, the fact that the observations slowly fall below the ^''Co line indicates 
that, as time goes by, more and more energy leaks out of the expanding remnant without contributing to the radiative 
output. Comparing the much more luminous SN 1998bw (magenta) to SN 1987A, one sees that the much higher ^^Ni 
production drives a more luminous peak, but the lower ratio of kinetic energy to ejected mass from the stripped progenitor 
results in an inefhcient trapping of the radioactive energy released, and the observed light curves falls rapidly below the 
energy release rate. Even more luminous SN 2006gy (magenta stars) would require > 7Mq of ^^Ni to reach the observed 



peak. However, as can be seen from comparison between the expected slow decay and the deep non-detections at late times 
(magenta inverted triangles), the observations are inconsistent with a radioactively driven evolution. Instead, this event 
is now understood as resulting from strong CSM interaction [7] . f8] . [14] . l43l . I44| . l45l . though of an unconventional sort with 
the interaction region initially opaque and invisible [57|, and the late-time luminosity came from a reflected-light ccho[8]. 
Finally, as we have reported here, the evolution of SN 2007bi is fully consistent with a PISN of an extremely massive star, 
producing several solar masses of ^®Ni (and enough ejecta to trap the radioactive decay energy), driving the light curve out 
to very late times, in perfect concordance with theoretical ^®Co decay. 



(5) Spectroscopic modelling 



Photospheric spectral fitting: 

In Fig. 1 we present an automatically derived [55] SYNOW|59j fit to our data. The 
best automatic fit derived agrees with our best manually derived attempts, and indi- 
cates a photospheric velocity fph = 12, 000 km s^^. Prominent lines of iron, calcium, and 
magnesium are seen, and flux depression in the blue side of the spectrum results from 
iron-group element (Fe, Co, Ni) line blends, as typically found for Type I SNe. Lines 
of neutral oxygen and sodium, which are often prominent in early Type Ic SN spectra, 
appear remarkably weak. Independent analysis by ref. [19] arrives at very similar results. 

Nebular spectral fitting: 

In Fig. 2 (bottom) we compare our Keck nebular spectrum with models derived 
using a well-tested nebular modelling code[17], operated in a single- zone mode. The code 
calculates radioactive excitation and nebular emission cooling using extensive line lists 
and constrains the amounts of radioactivity and the mass of the various elements. Since 
cooling effects at all wavelengths (extending outside of the optical window) are considered, 
elements lacking strong optical nebular lines are also constrained. Since the same code 
was used to study SN 1998bw[17], the relative results (SN 2007bi vs. SN 1998bw) are 
quite robust, and indicate that SN 2007bi produced a factor of ^ 10 more ^^Ni than SN 
1998bw. 



(6) The spatial structure of the ejecta at late times 

In Fig. 6 we present a schematic illustration of the apparent geometrical distribution 
of the ejecta. The main feature is that radioactive ^^Ni (decaying into ^^Co at late times) 
appears to be centrally concentrated, and not to have been mixed all the way out into the 
outermost layers of the envelope. It thus illuminates mainly the more slowly expanding, 
heavy-element-rich inner ejecta. The outermost, faster layers are not illuminated at late 
times, explaining the slower velocity of the material emitting the nebular spectra, and 
probably its composition, which appears to be depleted in C, O, and Mg relative to Fe. 
The outermost helium layer would probably lie even farther outside, well away from most 
of the ^^Ni synthesized, and thus does not contribute to either early (photospheric) or 
late (nebular) spectra. This simple spherical scheme appears to explain all available data, 
including the following: 

(a) The high mass estimated from light-curve modelling, which is sensitive to all of the 
material contributing to the opacity at the photospheric phase, including the outer, faster 
shells, compared to the lower total mass derived from nebular spectroscopic modelling, 
sensitive only to slower, inner shells which are highly enriched in radioactive material. 

(b) The composition derived from the nebular spectrum which, relative to iron, is 
depleted in lighter elements (C, O, Mg) that are more abundant in the outer layers of the 
envelope. 

(c) The lack of helium lines at all times, which are segregated from the energizing 
^^Ni and are thus not excited[41j. 

We note that the analysis of ref. [19] does not show evidence for asphericity in late- 
time spectra of SN 2007bi as seen in many other Type Ib/c events pO| ED E21 E3]- Along 
with the apparent suggestion that ^^Ni is not well mixed, the data probably argue against 
a bipolar /jet-driven explosion model as proposed for normal and GRB-related SNe Ib/c. 




Figure 6: 

A schematic illustration of the ejecta geometry. ^^Ni (decaying to ^^Co) is most abun- 
dant in the core, and emits positrons and gamma rays that excite the surrounding ma- 
terial. The outer layers of the envelope, which are expected to be dominated by lighter 
elements (C, O, Mg) and where any helium must reside, are not well mixed with the 
radioactive elements and thus remain unilluminated by hard radiation at late times, and 
do not contribute to the nebular spectrum (nor to the derived mass and composition from 
the analysis of these data). 
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